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ABSTRACT 

Context. L-type ultra-cool dwarfs and brown dwarfs have cloudy atmospheres that could host weather-like phenomena. The detection 
of photometric or spectral variability would provide insight into unresolved atmospheric heterogeneities, such as holes in a global 
cloud deck. Indeed, a number of ultra-cool dwarfs have indeed been reported to vary. Additional time-resolved spectral observations 
of brown dwarfs offer the opportunity to further constrain and characterize atmospheric variability. 

Aims. It has been proposed that growth of heterogeneities in the global cloud deck may account for the L- to T-type transition as 
brown dwarf photospheres evolve from cloudy to clear conditions. Such a mechanism is compatible with variability. We searched for 
variability in the spectra of five L6 to T6 brown dwarfs in order to test this hypothesis. 

Methods. We obtained spectroscopic time series using the near-infrared spectrographs ISAAC on VLT-ANTU, over 0.99-1.13 /mi, 
and SpeX on the Infrared Telescope Facility for two of our targets, in J, H and K bands. We search for statistically variable lines and 
correlation between those. 

Results. High spectral-frequency variations are seen in some objects, but these detections are marginal and need to be confirmed. 
We find no evidence for large amplitude variations in spectral morphology and we place firm upper limits of 2 to 3 % on broad-band 
variability, depending on the targets and wavelengths, on the time scale of a few hours. In contrast to the rest of the sample, the T2 
transition brown dwarf SDSS J1254-0122 shows numerous variable features, but a secure variability diagnosis would require further 
observations. 

Conclusions. Assuming that any variability arises from the rotation of patterns of large-scale clear and cloudy regions across the 
surface, we find that the typical physical scale of cloud cover disruption should be smaller than 5-8 % of the disk area for four of our 
targets, using simplistic heterogeneous atmospheric models. The possible variations seen in SDSS J1254-0122 are not strong enough 
to allow us to confirm the cloud breaking hypothesis. 

Key words. Stars: low mass, brown dwarfs - Stars: atmosphere - Technique: spectroscopic - Stars: individual: 
2MASS J08251968+21 15521, SDSS J125453. 90-012247.4, 2MASS J12255432-2739466AB, 2MASS J15344984-2952274AB, 
SDSS J162414.37+002915.6 
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1. Introduction 

L/T transition brown dwarfs are an informal class of sub-stellar 
objects comprising the latest L-type and the earliest T-type 
dwarfs (roughly L8 to T4). Their near-infrared spectra may ex- 
hibit both CO and CH4 absorption and they have J - K colours 
intermediate between those of the red late-type L dwarfs and 
the blue mid-T dwarfs (e.g. < J - K < 1.2). The spectra of 
L dwarfs are characterized by absorption due to metal hydrides 
(FeH, CrH) and strong alkali (K, Na, Rb, Cs) lines, absorption 
bands of water and CO, and the grey opacity of a cloud deck 
(Kirkpatrick et al. 119991 Martin et al. I19991 ). These clouds are 
made of condensates of refractory elements such as Al, Ca, Fe, 
Mg, Si, as well as Ti and V whose oxides dominate the opti- 
cal spectral morphology of M dwarfs, but which are depleted 
from cooler atmospheres as they condense into dust grains. The 
vertical distribution of the clouds is a complex issue, but a bal- 
ance between sedimentation (the equivalent of rain in the case 
of water) and upward mixing likely controls their vertical profile 
(e.g., Ackerman & Marley 2001). As the atmosphere becomes 
cooler in late-type L dwarfs, the grey clouds gradually settle be- 
low the photosphere and the emergent spectrum becomes more 
strongly affected by molecular band opacities. In the next cooler 
spectral class, T, CO completes its reduction to CH4 (Burgasser 
et al. |2002al Geballe et al. |2002| l. The latter produces strong IR 
absorption, along with water and molecular hydrogen (collision- 
induced absorption). Synthetic spectra predicted by atmosphere 
models with a discrete cloud deck reproduce well the colours 
(Ackerman & Marley 120011 Burgasser et al. 1200 2b). the near- 
infrared spectra (Cushing et al. l20051 l. and the mid-infrared spec- 
tra (Roellig et al. [2004) of the L dwarfs. Likewise models with 
no grain opacity fit well the spectra of mid- to late-type T dwarfs 
(Cushing etal. 12008b . 

However the rapidity of the transition between the cloudy 
L dwarf atmospheres and the clear T dwarf atmospheres is not 
well understood. Proposed mechanisms have included horizontal 
fragmentation of the global cloud deck (Burgasser et al. [2002b), 
a rapid increase in the cloud particle sedimentation efficiency 
(Knapp et al. 2004), or the sinking of a very thin cloud layer 
(Tsuji et al. 2004). The case for a dynamic mechanism, as 
opposed to a gradual sinking of a uniform cloud deck, was 
strengthened by recent measurements of geometric parallaxes 
for a significant number of brown dwarfs (Vrba et al. [2004), 
and the observations of brown dwarfs over a larger wave- 
length range of their spectra, particularly in the mid-infrared 
(Golimowski et al. 2004). Effective temperatures derived in- 
dependently of atmospheric modeling assumptions reveal that 
the L/T transition takes place over a small effective tempera- 
ture range (about 200 K), as was qualitatively first noticed by 
Kirkpatrick et al. ( |1999l >. This is difficult to reconcile with a 
gradual settling of the clouds. It is also marked by a brightening 
in the J band with increasing spectral type (Dahn et al. 2002; 
Knapp et al. l2004l their Fig. 8 and 9). Burgasser et al. (I2002bb 
interpreted this rapid brightening, and the corresponding J - K 
colour change, as an indication of the fragmentation of the global 
cloud cover in the L/T transition brown dwarfs' atmospheres. 

One possible way to test the horizontal fragmentation hy- 
pothesis would be to search for variability over a wide range of 
spectral types. Large-scale, patchy clouds, like those on Jupiter, 
would be expected to create photometric variations (Gelino & 
Marley 2000) as the structures rotate across the disk. If the sig- 
nature of such phenomena is preferentially observed in the L/T 
transition objects this would strengthen the case for cloud frag- 
mentation. If the typical length scale of horizontal cloud patches 



are small, however, little to no photometric or spectral variation 
would be expected. This would also be the case for stable stripes 
along the brow dwarf parallels. Thus detection of variability sup- 
ports the cloud fragmentation hypothesis, but non-detection can- 
not rule it out. Here we do not put forward any hypotheses on 
how such heterogeneities could arise. We aim simply to test if 
the surface is heterogeneous on large spatial scales. 

Other explanations for the L/T transition enigma have been 
suggested. Since we started this project, a number of L/T brown 
dwarfs have been resolved into binaries composed by a late-L or 
early-T primary, and a fainter T companion (e.g. Reid et al. 2001 
Burgasser et al. 2005 j Liu et al. 2006). Although the statistics are 
still not conclusive, the seemingly higher rate of binaries among 
L/T transition brown dwarfs compared to other brown dwarfs 
has prompted Burgasser et al. ( 2005 , 2006b ) and Liu et al. 2006 ) 
to speculate that binarity could explain the peculiar relation be- 
tween effective temperature and spectral energy distribution ob- 
served over the transition. In the case where most L/T transi- 
tion objects would be composed of earlier- and later-type com- 
ponents, we would not expect to find more variability among 
transition objects than among the populations from which the 
components are drawn, taking into account the binary effects. 
Similarly, our data do not allow to search confirmation of any 
mechanism that would not produce surface heterogeneities, in- 
cluding uniform and gradual sinking, or thinning (e.g., Knapp 
et al.[2004, Tsuji et al. [2004b of the clouds. 

There have been a number of reports of brown dwarfs ex- 
hibiting variability. Tinney & Tolley dl999l l, Bailer-Jones & 
Mund t Q99p |20TH1 |. B ailer-J ones & Lamm ( 120031 Clarke 
et al. d2002bl ) and Koen (2003 ) studied samples of early to mid 
L dwarfs in the / band and reported variability detections for 
a fraction of their targets. Gelino et al. (120021) conducted opti- 
cal and near-infrared photometry follow-up of a dozen of brown 
dwarfs, and found several variations in their targets' light curves. 
Enoch et al. (2003) conducted a A', -only monitoring program of 
nine L and T dwarfs, including three L/T transition dwarfs. At 
the 2-cr level, they found three variable dwarfs, including one 
Tl dwarf with a possible 2.96-h period, at the 89 % confidence 
level (C.L.), and did not find a correlation between variability 
and spectral type. Their average detection limit of sinusoidal 
variations is 12% (99% C.L.). Clarke et al. d2002al l reported 
that Kelu-1 (L3, Geballe et al. 120021 ) varies photometrically in 
their 2000 optical observations, with a 1.8 h period, but this was 
not observed again in 2002 (Clarke et al. 120031 ). Bailer- Jones & 
Lamm (2003 ) studied the effect of second-order extinction resid- 
uals in photometric variability observations. These residuals are 
due to the variable, wavelength-dependent telluric absorption of 
H2O, CO2 and O3 and to the very different colours of the brown 
dwarfs and the comparison stars. They estimate these residuals 
to be of the percent level, and cannot generally be removed. Our 
spectroscopic observations are mostly not affected by this effect. 
Only deep, non-resolved lines in both the target spectrum and 
the telluric absorption spectrum could be affected. 

Ha has been monitored in two L3 dwarfs (Hall 2002 Clarke 
et al. 120031 and an L5 dwarf (Liebert et al. 120031 , and the 
emission was found variable in all three objects. Burgasser 
et al. (12002c) studied the peculiar, Ha emitting T6.5 brown 
dwarf 2MASS J12373919+6526148 using optical spectroscopy 
and /-band photometry. A tentative shift in the Ha line was de- 
tected but no photometric variability, at the =s 2.5 % level. A se- 
ries of eleven 20-min spectra of the L8 dwarf Gl 584C revealed 
"subtle" qualitative variations in the 9250-9600 A water absorp- 
tion band, although a telluric origin cannot be excluded, and red- 
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ward of the CrH 8600 A bandheaad (Kirkpatrick et al. |2001] ). 
Near-infrared spectroscopic data are more scarce. Nakajima 
et al. (120001 1 cautiously reported a variation in the water band 
at 1.53-1. 58 pm in the spectra of the T6 SDSS J 1624+0029, us- 
ing the Subaru telescope, although over only 80min of time. 
The work most similar to this has recently been published by 
Bailer- Jones (2008 ). Two out of four ultracool dwarfs have been 
found to show correlated spectral numerous variations with 13- 
nm resolution, often (but not always) associated with one or 
more known features, particularly water and methane. 

Finally, Morales-Calderon et al. ((2006 ) searched for photo- 
metric variability in three late L dwarfs (DENIS-P J0255-4700, 
2MASS J0908+5032, and 2MASS J2244+2043) with the IRAC 
photometer on Spitzer Space Telescope. Two out of the three 
objects studied exhibited some slight modulation in their light 
curves at 4.5, but not 8.8 pm. 

In short, despite numerous variability searches with a great 
variety of instruments and technique, there has yet been no 
definitive detection of luminosity variability in a field brown 
dwarf. Here we report on a search for spectroscopic variability 
in a selection of L/T transition objects by the Clouds collabora- 
tion (Continuous Longitude Observations of Ultra-cool Dwarfs, 
Goldman et al. 2003 Q In the following we present our obser- 
vations (Section|2]i, and their reduction and how we search for 
variability (Section[3]). In Section|4]we detail for each target the 
variability information we obtained from our time series, and 
present variability upper-limits. Finally in Section|6]we interpret 
our results and link them to atmospheric physical processes. 

2. Observations 

2.1. Target selection 

We chose our targets according to several criteria. The targets 
had to sample the L/T transition spectral types, from late L to 
early T. The J band magnitude of the target had to be brighter 
than 15.5, in order to obtain a high signal-to-noise ratio (SNR) 
over most of the wavelength range, with a good time sampling. 
SDSS J1254-0122 was a natural target as we had detected pho- 
tometric variations in J, H and MKO K bands in April 2002 
(Goldman et al.|2003]). 

We selected brown dwarfs having a nearby star within 90 " . 
VLT/ Isaac has a long slit of 2' that allowed us to simultaneously 
observe a second star (hereafter: comparison star). This simpli- 
fies and strengthens the analysis as the comparison star can be 
used to remove most of the sky absorption variability that af- 
fects the target spectrum. We required the comparison star to be 
brighter than the target in the J band. The actual ratios of the 
comparison star flux to the target flux at 1 .08 pm vary between 2 
(SDSS J1254-0122) and over 30 (2MASS J0825+21 15). These 
constraints limit the number of possible targets to just a few. 

A T6 brown dwarf, SDSS J1624+0029, was selected to com- 
pare L/T transition object variability characteristics with later- 
type objects, and because it was reported variable based on near- 
infrared spectroscopic observations (Nakajima et al. 2000). We 
also observed the mid-T dwarfs 2MASS J1534-2952AB (T5) 
and 2MASSJ1225-2739AB (T6), when observing conditions 
prevented us from monitoring other targets. These two binaries 
are not resolved in our observations. 

Table Q] presents the final list of targets with a short de- 
scription, including binarity and previous variability informa- 
tion. In this article, we use near-infrared classification schemes 

1 http : //astronomy . nmsu . edu/CLOUDS/ 



of Burgasser et al. (2006a) and of Geballe et al. (120021 ) in the 
case of 2MASS J0825+21 15. The target names are truncated to 
the first four digits of the right ascension and declination. 

2.2. VLT/ISAAC spectra 

We used Isaac, the Infrared Spectrometer And Array Camera 
mounted on VLT-Antu, during the nights of 2003, April 
16 and 17. Table|2] shows the observing log. The low- 
resolution grism covers the wavelength range of 0.99-1.1 3 pm 
with a resolving power of R = ■A « 550. Both nights 
were partly cloudy. Relative extinction and slit losses for 
each target were up to 0.55 and 1.0 mag for April 16 and 
17, respectively. When we observed the Northern targets 
(2MASS J0825+21 15, SDSS J1254-0122, SDSS J1624+0029), 
the winds, from the North, made the seeing highly variable. We 
observed 2MASS J1534-2952AB and 2MASS J1225-2739AB 
when stronger winds prevented us from observing towards the 
North. Most observations were conducted using the l"-wide slit 
(corresponding to a 7-pixel resolution element), except at the 
beginning and the end of the second night when we used the 
l'.'5-wide slit (10-pixel). The site seeing during the first night 
varied between 076 to 172, increasing sharply at the end of 
the SDSS J1254-0122 observations and the beginning of the 
2MASS J1225-2739AB observations. During the second night, 
the seeing increased slowly from 0795 to 172. 

Spectra were chopped between two positions, A and B, sep- 
arated by 20". In order to minimize the impact of bad pixels and 
the systematic effects of inaccurate flat fielding, the A and B po- 
sitions were randomly moved around each center by up to 5". 
We chose a single integration time of 2.5 min to minimize sky 
brightness fluctuations between the A and B exposures. All tar- 
gets were observed with a comparison star in the slit, and this 
constraint sets the slit angle. The parallactic and slit angles are 
indicated in Table|2] The centering of the stars was typically 
checked after two hours of observations. The centering proce- 
dure causes the grism wheel to be moved; no other changes in 
instrument configuration occurred during each observing night. 

After the last scientific observations of each target, except 
for 2MASS J0825+21 15 at the beginning of the second night 
and 2MASS J1534-2952AB at the end of the second night, we 
took a lamp flat, a bias, and a series of arcs for wavelength cali- 
bration. A spectrometric standard, the AO V star HD 1 1 174, was 
occasionally observed (once the first night, twice on the second 
night). As we are not interested in absolute spectra, we only took 
these for future reference. The stacked VLT spectra of April 
16 (April 17 for SDSS J1624+0029), wavelength- and flux- 
calibrated, are shown in Fig.Q](top). The SNR of the spectra over 
1.0-1.1 pm ranges from 100 to 200 (75 for SDSS J1624+0029). 
The features seen over 1 to 1.1 pm are mostly real. Several 
features in the SDSS J1254-0122 spectrum match those of the 
NIRSPEC spectrum published by McLean et al. d2"0"03"i Others 
are barely significant and may be due to systematic errors or 
noise. 

2.3. IRTF/SpeX spectra 

We obtained spectra of 2MASS J0825+ 21 15 and 
SDSSJ1254-0122 using SpeX (Rayner et al. 120031 the 
medium-resolution cross-dispersed spectrograph mounted 
on the 3.0 m NASA Infrared Telescope Facility (IRTF), 
on 2003 April 18 to 21 (UT). Tabled shows the ob- 
serving log. The stacked IRTF spectra of April 18 and 
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Table 1. Summary of targets and references. 



Target 


SpeX Sp.T. 1 


J 


Mi 


separation 2 


discovery 


variability 



2MASS J08251968+21 15521 


yes 


L6 


14.89 


14.84 = 


t 0.08 3 




4 


Kirkpatrick et al. (12000b 




2MASS J12255432-2739466A 




T6 


15.50 


14.92 = 


t0.14 3 


282 


± 5 mas 


Burgasser et al. (19991 




2MASS J12255432-2739466B 




T8 


16.85 


16.27 = 


t0.14 3 






Burgasser et al. ( 2003 1 




SDSSJ125453.90-012247.4 


yes 


T2 


14.66 


14.00 = 


t 0.06 5 




6 


Leggett et al. d2000> 


J, H, K 8 


2MASS J15344984-2952274AB 




T5.5 


14.60 5 


13.94 = 


t 0.06 7 


65 : 


± 7 mas 


Chiu et al. U0061 




SDSS J162414.37+002915.6 




T6 


15.41 3 


15.13 


+0.06 3 
-0.05 






Strauss et al. l|1999|> 


H 2 0" 



All targets were observed with Isaac 

1 All spectral types from Burgasser et al. d2006at except for 2MASS J0825+21 15 (Geballe et al. l2"0"02l >. 

The photometric type of 2MASS J12255432-2739466B was determined from photometry, and is hence uncertain. 

2 Binary separations, 2MASS J magnitudes and spectral types for each component from Burgasser et al. (2003 ). 

3 2MASS J photometry f rom V rba et al. J2004t . 

4 Unresolved (Bouv et al. l2003l . 

5 MKO J magnitude from Kn app et a l. J2004t . 

6 Unresolved (Burgasser et al. 2006b, Goldman et al., in prep). 
1 MKO J magnit ude fro m Tinney et al. ( [20031 . 

8 Goldman et al. (2003 ) reported photometric variations in J, H and K bands. See Section [5T2l 

9 Nakajima et al. (2000) reported spectroscopic variations in the water bands over 1.53-1.58 //m. 



Table 2. Isaac observing log. The reported seeing values are those of the DIMM seeing monitor; seeing measured on the frames 
show very similar values. 









slit 






# of spectra 


slit 


parallactic 


Date 


Target 


Start UT 


width 


site seeing 


Duration 


[# taken] 


angle (°) 


angle (°) 


April 16 


2MASSJ0825+2115 


00:10:49 


1" 


a' 82-0".' 55 


70min 


20 


114 


167-148 




HD 11174 


01:30:20 


1" 


or 55 


1 min 


2 


179 


144-49 




SDSSJ1254-0122 


01:42:05 


1" 


0780-1705 


208 min 


66 


141 


45-145 




2MASSJ1225-2739AB 


05:30:33 


1" 


(y.'80-l'.'25 


90 min 


30 


54 


90-100 




2MASSJ1534-2952AB 


07:49:22 


1" 


0795-1'.' 40 


144 min 


50 


98 


76-97 


April 17 


2MASSJ0825+2115 


00:17:01 


1.5" 


0780- 1'.' 10 


56 min 


20 


114 


164-149 




HD 11174 


01:23:52 


1.5" 


07 90-17 1 


15 min 


8 


179 


44-41 




SDSSJ1254-0122 


01:50:17 


1.5" 


07 80-17 1 8 


103 min 


36 


141 


50-12 




SDSSJ1254-0122 


03:49:06 


1" 


0785-1722 (1740) 


163 min 


44 [52] 


141 


182-124 




HD 11174 


06:40:53 


1" 


1702-1724 


2 min 


4 


179 


131 




SDSS J 1624+0029 


07:26:58 


1.5" 


0772-1740 


150 min 


30 [42] 


30 


178-127 



19 for 2MASSJ0825+2115 (April 18, 20 and 21 for 
SDSS J1254+0122), wavelength- and flux-calibrated, are 
shown in Fig.Q] (bottom). Observing conditions were non- 
photometric on all nights and observations were often halted 
due to thick cloud coverage. The seeing was typically ~078 in 
the //-band. We used the ShortXD mode with the 078-wide slit 
that covers from 0.8-2.4 /urn at R a 750 in a single exposure. 
All observations were conducted at the parallactic angle to 
minimize differential slit losses. No comparison stars were 
observed simultaneously because of the 15 "slit length. 



For each target, we obtained a series of 2-min exposures. 
During the series, the telescope was nodded between two po- 
sitions along the 15" slit to facilitate the subtraction of the 
dark current and sky background in the data reduction process. 
Observations of an AO V standard star were also acquired to cor- 
rect each target for telluric absorption. The airmass difference 
between the targets and the AO V stars was always less than 0. 1 . 
Finally, a set of calibration frames including flat field and argon 
arc exposures were obtained for each target-standard star pair. 



Table 3. SpeX observing log. All slit angles were parallactic. Slit 
width was 078. Spectroscopic standard observations, between 
sets of target observations, are not reported. 



Date 


Target 


Start UT 


Duration 


spectra 


April 17 


2MASSJ0825+2115 


06:19:38 


42 min 


18 


April 18 


2MASSJ0825+2115 


05:35:48 


18 min 


8 




SDSSJ1254-0122 


08:32:05 


42 min 


18 






09:31:41 


42 min 


18 






10:25:49 


28 min 


9 






12:30:03 


18 min 


8 


April 19 


2MASSJ0825+2115 


05:33:27 


20 min 


10 






06:15:03 


21 min 


9 






07:43:54 


42 min 


18 


April 20 


SDSSJ1254-0122 


10:07:23 


21 min 


10 






10:45:46 


42 min 


20 






10:50:29 


26 min 


11 






12:27:47 


22 min 


10 


April 21 


SDSSJ1254-0122 


10:43:02 


20 min 


10 






11:21:41 


41 min 


20 






12:19:58 


22 min 


10 
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3. Reduction 

3.1. VLT/Isaac spectra 

We reduce the spectra following the Isaac Data Reduction 
Guide 1.5. We remove the electric ghosts generated by the 
HAWAII detector using the ESO eclipse ghost recipe. The next 
steps are performed using a modified version of reduction pack- 
age Spextool v3.1 (Cushing et al. 2004). We subtract the A and 
B images to remove the sky background and bias, and the re- 
sulting images are flat-fielded. The creation of the flat field is a 
crucial step and its accuracy will be studied in Section [4~T| The 
residual sky background is subtracted, and the spectra are ex- 
tracted, using an aperture independent of wavelength. For each 
target, we wavelength calibrate each spectrum relative to a refer- 
ence sky spectrum, before we calculate the wavelength solution 
for this reference spectrum. We extract a sky spectrum corre- 
sponding to the A position from the B image, and shift the A 
object spectrum by a constant, using the OH lines of the sky 
spectra. Because this constant is usually not an integer, in units 
of wavelength-pixels, the flux values are linearly interpolated for 
the resulting fractional-pixel shifts. We estimate the errors intro- 
duced by this process to be about 1 % in rapidly changing parts 
of the spectrum (water bands) and much less elsewhere. For the 
two targets observed both nights, a shift of about 15 pixels oc- 
curred between the two nights of observation, much larger than 
shifts within a night. 

The reference sky spectra are calibrated to actual wave- 
lengths using the arc spectra taken right after the observations, 
and checked against OH sky lines, with a 3—4 A agreement. 
Absolute calibration is only required for line identification, not 
variability detection. 

The target spectra are then divided by the comparison star 
spectrum. Because we observed over 0. 99-1.1 3 //m where the 
Earth atmospheric transmission is good, the main effect of this 
correction is to remove terrestrial grey atmospheric absorption 
(clouds), instrumental response, and slit losses, mostly due to 
seeing variations. The comparison star, however, is useful to 
remove the water-band variations redward of l.lOjt/m. The re- 
sulting spectra are co-added by groups of six to eight (15 to 
20 minutes), depending on the target. It proved unnecessary to 
apply a median filter or clip the data, as most of the cosmic rays 
and bad pixels are removed by Spextool during the spectrum 
extraction. We choose to co-add six to eight spectra into stacked 
spectra, in order to get sufficient SNR, while maintaining sensi- 
tivity to short term variations. We search for variability on these 
15- to 20-min spectra, divided by the rest of the night data set. 
Finally, the overall flux level (independent of wavelength) may 
vary differently for the target and the comparison star, whenever 
the centering of both objects is not identical, so that the amount 
of flux losses differs. Depending on targets (and observing con- 
ditions), we sometimes have to normalise each spectrum over a 
broad wavelength region, when the stability obtained after cor- 
recting with the comparison star's spectrum is not good enough 
for our purposes. For homogeneity, we apply this normalisation 
to all VLT targets. To determine that region, we try to maximise 
its extend, in order to be sensitive to variations affecting large 
parts of our wavelength coverage, and to increase the accuracy, 
while avoiding zones of artefacts (ghosts) and low signal (water 
bands, order edges). The two latter requirements vary with tar- 
gets, and the exact values are given in the legend of Fig.|2j|5]|2l[8] 
and|9] The resulting regions are 38-nm wide (27 % or our wave- 
length coverage, for SDSS J1624+0029) or wider. Searching for 
broad band photometric variations (outside atmospheric water 



bands) is in any case better performed with imaging. The disper- 
sion of the stacked spectrum flux values is used as a measure of 
the (total) error (see 14.21 ). The spectra presented in Fig. [2] [HE] [8] 
and|9]are smoothed over 1 1 pixels, or 1.5 times the actual reso- 
lution obtained with the l"-wide slit. 

3.2. /RTF/SpeX spectra 

The data are reduced using Spextool (Cushing et al. 20043, the 
IDL-based data reduction package for SpeX. Pairs of exposures 
taken at the two different slit positions are first differenced to 
remove the dark current and sky background (to first order) 
and then flat fielded. The spectra are extracted, using an aper- 
ture independent of wavelength, and wavelength calibrated after 
subtracting any residual sky background. The raw spectra are 
then corrected for telluric absorption using the AO V star spectra 
and the technique described by Vacca et al. (2003). Finally, the 
telluric-corrected spectra from the different orders are merged to 
create continuous ~ 1-2.4 fim spectra for each target. 

As for the VLT spectra, the dispersion of the flux values is 
used as a measure of the error (see !4.2l for a brief discussion of 
the drawbacks of this method). We also normalise each spec- 
trum over J (more exactly 1.19-1.32 /mi), H (1.44-1.75 /im) and 
K (2.02-2. 19 fim), and perform our variability search indepen- 
dently in each region of the spectrum. 

The spectra presented in Fig. [3] and [4] are smoothed over 
13 pixels. 

3.3. Slit losses, refraction and the like 

Varying observing and instrumental conditions (such as move- 
ments of the slit relatively to the target, or seeing variations) 
occurring during the observations could be a problem for our 
study. Because we aim at a great precision (better than 1 %) of 
our relative spectroscopy, we need to pay special attention to 
flux losses that could mimic intrinsic variations. The effects de- 
scribed in this section vary smoothly with wavelength. They will 
not affect the analysis of high-frequency variations, but mostly 
the calculations based on the spectral indices of Section lrxTl 

Indeed, the slit cuts the wings of the point-spread function 
in variable amounts depending on the wavelength. Seeing varies 
with wavelength, according to: s oc A~ 2 . Over the Isaac wave- 
length range, this means that the flux loss is smaller by 1 % 
at 1.13/vm relative to 0.99 //m (assuming 1" seeing and 1" slit 
width, and the stars to be centered on the slit). This effect, how- 
ever, affects both the target and the reference star equally and is 
therefore removed. Within the SpeX J, H and K s bands, the ef- 
fect is 2 % to 1 %, for 0.8" to 1.2" seeing (0.8" slit width). Most 
of this is removed using the reference stars, which were observed 
every 20 to 40 min, hence under similar seeing conditions. 

The stellar centroid also shifts with wavelength due to atmo- 
spheric refraction, depending on the hour angle. For the IRTF 
data, we observed at the parallactic angle, so that the star was 
equally centered on the slit at all wavelengths. The main limi- 
tation here is the wavelength-dependent variations of the atmo- 
spheric transparency between two observations of the reference 
star. 

For the VLT data, the slit angle was set by the nearby refer- 
ence star (see Table|2](. We are now concerned with differential, 
wavelength-dependent flux losses between the target and the ref- 
erence star. Slit translation movements and seeing variations, if 
the target and the reference star are not equally centred in the slit, 
and slit rotation, will result in wavelength-dependent flux losses, 
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due the wavelength-dependent centroids. (If the two objects are 
equally centred on the slit, flux losses will be similar and can- 
cel out.) Slit rotations are not expected with Isaac Occasional 
flexures are not known to produce significant effect with our slit 
widths. For two targets, we performed a second acquisition to 
check orientation after a long observing sequence, and found 
no significant rotation. The VLT unit telescopes are altazimu- 
tal but none of our targets were observed closer than 17 deg to 
the zenith, while 10 deg is considered a safe limit. 

We nevertheless investigate the possibility of a colour effect 
of a small rotation. We estimate the differential refraction shift 
using the procedure of Marchettfl The maximal relative shift 
predicted between 0.99 /mi and 1.12 /jm is 0.05 ", when compar- 
ing observations taken at the minimal and maximal airmasses of 
each target. 

In practice, our spectral index svlt (see 16. Q does not show 
any correlation with the normalisation factor, for any target. We 
conclude that refraction and slit losses do not add systematic 
effects to our measurement. They may increase the statistical 
noise, although this will only be perceptible as we add fluxes 
of many pixels. 

4. Analysis 

4.1. General considerations on the VLT/Isaac spectra 

The correction by the comparison star allows us to obtain a 1- 
3 % precision in most co-added spectra, before normalisation. 
This good result proves the usefulness of observing a compari- 
son star simultaneously with the target. The Isaac flat fields are 
known to be very stable over several nights for high frequency 
scales. The set up (grism, etc) wasn't changed during the ob- 
servations. Only the grism wheel was rotated to center the star in 
the slit about every second hour. The instrument is on a Nasmyth 
focus which reduces the distortions, such as wavelength shifts. 
Flat field errors are an important source of systematic errors. We 
compare the flat fields taken during the same night, with a few 
hour gap. A slope of about 0.5 % and broad variations of similar 
amplitude are observed. We do not remove these effects. They 
are always insignificant compared to the photon noise of the 
brown dwarfs. When we stack the spectra of the brightest com- 
parison stars, we find that the flux variations over 0.99 to 1 . 1 1 fim 
are 0.43%, 0.52% and 0.45% (per pixel) for the compari- 
son stars of 2MASSJ0825+2115, 2MASS J1225-2739AB and 
2MASS J1534-2952AB, respectively, with no slope larger than 
0.05% per 0.1 jum. 

Flat fields taken the same night, but for different targets or 
different slit widths, do show a 5 % variation at 1 .09 fim. As this 
is seen over the whole chip, this effect is mostly removed by di- 
viding the spectra by a median spectrum. Because we generally 
used flat fields taken right after the observations, it is possible 
that this effect does not affect our data. Consequently, system- 
atic effects are likely to be less than, or occasionally equal, to 
0.5%. 

The flat fields between two nights vary much more, due to 
day-time instrument configuration changes for routine calibra- 
tion (grism position, slit, flat field lamp temperature). Although 
some broad-band effects can be removed by fitting a 3-5-degree 
polynomial, the systematics increase significantly, and we did 
not correct for these variations. 

A "ghost" of each spectrum, scaled by about 1 .5 %, is present 
at ~1024 - y pixels where y is the position of the spectrum along 

2 eso.org/gen-fac/pubs/astclim/lasilla/diffrefr .html 



the spatial axis, over 525-635 rows in the dispersion axis (about 
1. 055-1. 086yt/m). This is probably an internal, optical reflection, 
or an electric cross-talk. The comparison star's ghost overlaps 
with a spectrum (comparison star or target) for about 19 spectra, 
which are only used over the unaffected region of the spectrum. 
About 20 ghosts overlap with a spectrum of the correspond- 
ing dithered frame. We examined those spectra and masked the 
affected wavelengths if the distance between the ghost of the 
brightest comparison star and the target spectrum is smaller than 
2". We could have subtracted a different image, but that would 
have increased the time difference between the pair of images 
and the contribution to the noise of the sky subtraction. The ef- 
fect of the target ghost on the faintest comparison star is 0.7 % 
at most and is neglected. 

4.2. Flux uncertainties 

In order to estimate the significance of possible variations, we 
need to know the noise associated with our measurements. We 
use the dispersion observed over 15 to 20 minutes (Isaac) and 
over 16 to 40 minutes (SpeX). This dispersion may include short- 
term intrinsic variability. It also includes any remaining system- 
atics affecting these short time scales. We would not be sensi- 
tive to 1 % variations over a few minutes, due to the limited 
SNR, anyway. Given the small number of spectra from which 
the dispersions are measured (six to eight for Isaac, eight to 20 
for SpeX), the uncertainties on the dispersion are large. The un- 
certainties are however reduced when we smooth the spectra in 
wavelength. 

4.3. Method 

We use coherent methods to detect variations and to calculate 
variability upper limits, in the case of non detections. The meth- 
ods are coherent because we perform identical variability detec- 
tions on the observed spectra and, in order to determine the upper 
limits, on simulated spectra. We base our search on the^- 2 prob- 
ability distributions of the time series. For each wavelength (for 
each pixel of a reference spectrum), we fit the time-dependent 
fluxes by a constant flux and obtain a^ 2 . Using a Savitsky-Golay 
filter, the stacked spectra are previously smoothed to match the 
resolution element (resolving power of R — 500). We also look 
for lower frequency (broader) variations, and perform a sim- 
ilar search after smoothing the spectra to a resolving power 
of R = 100. We decide that a variation is significant when the 
X 2 value is compatible with the random distribution hypothesis 
with a probability lower than 1 % (99 % of confidence level). For 
each target, we comment on those variations in Section[5] 

We calculate the integrated flux around those features which 
appear to vary. We obtain a narrow-band photometric time se- 
ries, and search for linear correlation between the fluxes obtained 
for different features. We report in Table|5]those features which 
show correlated variations, with a probability higher that 99 % 
of confidence level. Because we perform these correlation mea- 
surements only on a limited set of lines, those that show varia- 
tions, we expect less than one spurious detection. Correlation in 
features widely separated in the spectrum (hence on the array) 
strongly argue that the variations are not of instrumental origin, 
but originate in either the target or the Earth's atmosphere. 

To calculate the variability upper limits, we proceed in the 
manner of Bailer-Jones & Mundt ( 12001 1 ). We simulate mock 
spectroscopic time series, adding noise according to the ob- 
served flux errors. For each wavelength, we add to each stacked 
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spectrum time-dependent variations with increasing amplitudes, 
until the variations are detected, under the procedure used for 
the actual data. We model the time-dependent variations with 
a sinusoid, positive function: 5F = | sin(wf + 0)|. This would 
be observed (having neglected the limb-darkening effect) in the 
on-average-favourable case of two anti-podean spots. Hence we 
assume that a spot is observable at all times, with a variable pro- 
jected area. This means that our upper limits are only valid while 
the target was being observed, not on average over the observ- 
ing run. This calculation is more appropriate for our limited time 
coverage. An hypothesis of numerous spots, simulated by a sum 
of randomly-phased sinusoids, would also reduce the expected 
peak-to-peak amplitude. 

Then, the smallest amplitude that triggers the variability de- 
tection is recorded as the variability upper limit for that wave- 
length. The mock variations are time-sinusoid, with random pe- 
riods (between 1 and 16 hours) and phases. We average these 
maximum amplitudes for various periods and phases. 

We perform those calculations for each wavelength inde- 
pendently. However, real variations would certainly be corre- 
lated in wavelength, because of the non-zero slit width and 
of the physical and chemical processes in the BD atmosphere 
causing the variations. As this latter correlation is largely un- 
known, our mock variations are not correlated with wavelengths. 
However, we stress that the upper-limits we calculate are fluxes 
of smoothed spectra; hence they refer to variations affecting the 
whole resolution element (for either of the 100 or 500 resolving 
powers). 

Finally, we introduce in Section |6TTl two spectral indices 
which maximize our predicted sensitivity to variations caused 
by changes in the cloud pattern, one for the VLT wavelength re- 

S 10n: S = 1.06-1.08^ ° ne fOT the IRTF data: S = 1.58-1.62^ and 

use them to obtain variability upper-limits. Although the index 
definitions were motivated by particular models, the variability 
upper limits arising from the measurements themselves are of 
course independent of the models. 

5. Results 

We now discuss each individual target, starting with earlier spec- 
tral types, according to the recipes described in Section l4~3l The 
spectroscopic time series for each target are presented in Figf2] 
to [9] In this section we report and comment on the detected sig- 
nificant variations. We discuss the upper limits we obtain in the 
next section. 

5. 1 . 2MASS J0825+21 1 5 (L6) 

On both April 16 and 17, we obtained a series of twenty 2.5-min 
Isaac spectra (see Fig.[2j. We combined them by groups of six or 
eight. The standard deviations of the 15-min spectra is 1.5 % per 
resolution element with peaks at 3 % in the strongest sky lines. 

At the resolution of R — 500, we find five features with sig- 
nificant variations, with width of one to three pixels, all dur- 
ing the second night (see Fig|2j. At R — 100, we find broader 
and more significant features: at 0.996 /an during April 16, and 
more clearly over 1.064-1. 068 //m, at 1.089 /mi and 1.1 11 /an 
(see Table|5]l. However, none of those features vary with any sig- 
nificant correlation with any other line. 

We also obtained 45 SpeX 2-min spectra on April 18 and 19, 
which we combined by groups of eight to 18, depending on the 
observing conditions (see Fig. [3]). The relative precision is about 
3-6 % in J band, 2-5 % in the H band and better than 2 % in the 



K band. We see one 1-nm wide variable line, at 1.140 /mi. The 
telluric water band absorption makes there a variability claim 
difficult, in the absence of simultaneous observations of a com- 
parison star. 

The Isaac 0.996/im feature could be related to the CrH 
0.99685 /mi bandhead (Kirkpatrick et al. [19991 and 0.997 /an 
FeH Wing-Ford absorption feature (Wing et al. 119771 Cushing 
et al. 120051) . In cool regions of the photosphere of late-L and 
T brown dwarfs, Fe is depleted from the gas to form clouds, 
that, along with the silicate clouds, have optical depths greater 
than unity at these wavelengths. However, when a hole in the 
cloud deck appears, deeper, warmer regions of the atmosphere 
are unveiled, so that the FeH absorption features might be ex- 
pected to re-appear and hence vary when the cloud deck evolves 
(Burgasser et al. 12002b j). Variations in the CrH and FeH band- 
head regions redward of 8600 A have previously been noticed 
in the slightly later L8 dwarf Gl 584C over a 3-month baseline 
(Kirkpatrick et al. 2001). However, in the absence of quantita- 
tive analysis of the significance and amplitude of the variations, 
it is difficult to compare our results. The 1 .089 //m and 1.111 /mi 
features are on the edge of water and methane absorption bands, 
but no other features appear to vary redward of those features, 
where absorption is larger, and our sensitivity unchanged (see 
Fig.[[3). 

5.2. SDSS J1 254-01 22 (T2) 

We obtained a series of 66 and 81 2.5-min Isaac spectra, on 
April 16 and 17 respectively (see Fig.[5j. We combined them 
by groups of eight, except the last ten of April 16 stacked in one 
single spectrum. The precision ranges from 0.8 % to 2 % at the 
strongest sky lines and on the order edges. 

At the resolution of R = 500, we find a number of features 
with significant variations, with width of one to four pixels. The 
most significant one is at 0.997 /mi, mostly on the first night, and 
less significantly on the second night. The variations at 1.046 /an 
and 1. 11-1.1 3 /an are well detected on both night, while other 
variations are only observed during one night. During the sec- 
ond night a total of 16 features appear to vary. None of those are 
strong features with deep unresolved absorption lines, so that tel- 
luric differential absorption between the target and the reference 
star could not explain them. At R — 100, we find no variable 
features during the first night but four features appear to vary 
on April 17, particularly at 1.033//m and around 1.104//m (see 
Table|5]l. Among those features, we find five pairs which are cor- 
related (see Fig|6] top). 

The most significantly variable features (see Table|5j are 
located close to but outside strong and variable telluric emis- 
sion lines, such as 0.997 /an (close to 0.995 /an), 1.031 /an 
(1.029 /an), 1.046 /an (1.042 /an). Moreover, no variations are 
detected on the emission lines themselves. The broadest feature 
1.11-1.13/an overlaps with no emission lines. The observing 
conditions of SDSS J1254-0122 in terms of airmass and atmo- 
spheric parameters were similar to, possibly more stable than, 
those of the rest of the samplfl 

We also obtained 146 2-min SpeX spectra on the nights of 
April 18, 20 and 21, which we stacked by groups of eight to 20 
(see Fig.|4|). At the resolution of R = 500, we find many features 
which show variability. Among those, a dozen of pairs show 
significant and convincing correlations. Correlated pairs cover a 
wide range of features, with no obvious chemical connections- 
except possibly water bands, which affect a large part of the near 

3 See http : //archive . eso . org/asm/ ambient-server 
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Fig. 3. IRTF spectroscopic time series for 2MASS J0825+21 15. The 16- to 36-min spectra are relative to the average of April 19, 
and offset. UT is the mean time since April 16, 2003, UT=0:00. Measurements with errors larger than ~8 %, as well as water bands, 
are not displayed for clarity. No variability information can be confidently obtained in those regions. All spectra are normalised 
independently in J, H and K. The dashed grey line at 1.140/im marks the feature that we report to vary. 
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Fig. 4. IRTF spectroscopic time series for SDSS J1254-0122. The 16- to 40-min spectra are relative to the average of each night, 
and offset. UT is the mean time since April 16, 2003, UT=0:00. See legend of Fig.[3]for details. In contrast to 2MASS J0825+21 15, 
many lines exhibit statistically significant variations. 
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Fig. 5. VLT spectroscopic time series for SDSS J1254- 0122, for April 16 (left) and 17 (right), 2003. The 20-min spectra are relative 
to the average of each night, and offset for clarity. UT is the mean time since April 16, 2003, UT=0:00. The SDSS J1254-0122 
comparison star is only 2.5 to 5 times brighter than the brown dwarf, and its variations are not over-plotted for clarity. All spectra are 
normalised over 1. 002-1. 088/mi. The bottom left panel shows the medium signal to noise ratio of SDSS J1254-0122 as a function 
of wavelength. See legend of Fig.|2]for more details. 



infrared spectrum. We see some variations around 1.58/mi in 
the water bands, as Nakajima et al. (2000) reported in the spec- 
trum of SDSS J1624+0029, although we find significant varia- 
tions over a smaller wavelength range. Finally, some features 
(1.079 /mi and around 1.112 /mi) vary in both data sets. 

SDSS J 1254- 122 h as been reported to vary in April 2002 
by Goldman et al. d20031 l. in MKO J, H and K photometry, at the 
5 % level, with no detectable colour variation, over 4 hours. Our 
IRTF spectra cannot bring constraints on these variations. 



By chance, a second comparison star, closer to the main 
comparison star, is visible in the slit, although it is not well cen- 
tered. The relative variations of this second comparison star are 
smaller than for our target and displayed in Fig. [7] We only used 
the brighter, well centered, star to correct the fluxes of both the 
target and the second comparison star. 

At the resolution of R = 500, we find no features with sig- 
nificant variations. At R = 100, we find two broad features at 
1 .046 /mi and 1 .068 /mi, which show variations but no correla- 
tion. 



5.3. 2MASSJ1534-2952AB (T5) 

We obtained 46 2.5-min Isaac spectra of 
2MASSJ1534-2952AB, at the end of the night of April 
16 (see Fig. IT). The site seeing was highly variable and in- 
creased with time during our observations. The precision ranges 
from 1.2-2% over 1.01 to 1.10 /mi to 4% on the order edges 
and on the strong sky lines. 



5.4. 2MASSJ1225-2739AB (T6) 

On April 16, we obtained a series of thirty 2.5-min Isaac spectra 
(see Fig. [8]), which we combined by groups of six. The precision 
ranges from 1.5-2% over 1.03 to 1.10 /mi to 5 % at 1.1 2 /mi. At 
the resolution of R — 500, we find three features with significant 
variations, with width of one pixel only. At R = 100, we find one 
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broader and more interesting line at 1.1 04 /mi, over 1.5 nm, but 
due to one stacked spectrum only (see Tabled- None of those 
features vary with any significant correlation with any other line. 

5.5. SDSS J1 624+0029 (T6) 

We obtained 42 2.5-min Isaac spectra of SDSS J1624+0029, at 
the end of the night of April 17, which we stacked by groups 
of six (see Fig.|9]l. The precision ranges from 2% over 1.03 
to 1 .09 fim to 5 % on the order edges and on the strong sky lines. 
The last 12 spectra are discarded due to high extinction and slit 
losses (totalling 1 mag). At the resolution of R = 500, we find 
four features with significant variations, of width one or two pix- 
els, which show no correlation. AtR = 100, we find one broad 
line at 1 .033 /mi with width 5 nm. There is a FeH absorption line 
at 1 .03398 /mi (Cushing et al. 2003); for such a cool object how- 
ever FeH is not expected to play a major role, and we do not see 
variations in the stronger 0.99-/mi band head. 

Nakajima et al. (2000) reported a variation in the water band 
at 1.53-1.58 jum using the Subaru telescope, although over only 
80 min of time. It is difficult to estimate the amplitude and sig- 
nificance of the variability in their data. We do not see significant 
variations at the edge of the water band absorption at 1 . 1 //m, in 
our Isaac spectra. 

6. Discussion 

6.1. Models of cloud fragmentation 

One of the proposed explanations to the reported brown dwarf 
variability and the rapid disappearance of the cloud deck at the 
L/T transition is the fragmentation of the cloud deck (Ackerman 
& Marley [200T1 Burgasser et al. I2002bl) . This is only one of 
several hypothetical mechanisms to produce variability. We 
compare our results with crude theoretical predictions based on 
that hypothesis. We assume that the flux coming from a partly 
cloudy brown dwarf is the sum of the flux predicted by either the 
Cloudy and Clear or Settl and Cond models, weighted by the 
cloud and clear coverage. Here Cloudy and Settl are the ver- 
tically constrained cloudy atmosphere models by Ackerman & 
Marley ( I2001D and Allard et al. d2003l 2006, priv.com.), respec- 
tively. Clear and Cond are the respective cloud-free models, 
which have similar, but not identical, assumptions built into their 
construction. This simulation should only be taken as a guide, 
in the absence of dynamic models. In Fig.[TT]and[T0] we show 
the relative effect of the appearance of a 5 % cloud (replacing a 
cloud-clear region) on a brown dwarf. The set of temperatures 
and parameters we use are intended to be representative of our 
targets rather than the best fitted values for each of our five tar- 
gets. We estimate that differences between our predictions due 
to change in substellar parameters are smaller than the uncer- 
tainties of the predictions themselves. For each temperature, we 
adopt as a reference cloud coverage expected by the atmosphere 
models. Hence, a 1 100 K brown dwarf is well reproduced by 
clear model, so we consider the spectroscopic variations intro- 
duced by the appearance of a 5 % cloud on a clear brown dwarf. 
Similarly a 1 600 K brown dwarf is well described by a cloudy 
model, and we show the spectroscopic variations introduced by 
the disappearance of a 5 % cloud on a fully cloudy brown dwarf. 
For the 1400K brown dwarf, such as SDSS J1254-0122, we 
compare a 55 % covered atmosphere to a 50% one. How does 
the actual ratio of cloudy to clear areas change the expectations? 
The importance of using different reference cloud coverages are 
illustrated in FigQl] A similar study over 1-2.5 /mi shows the 



same trend, except in the water bands (1.4/zm, 1.9 /mi) where 
we have little or no signal. The relatively larger variations seen 
in the water band of cooler brown dwarfs are due to the smaller 
flux escaping at those wavelengths. 

In the VLT 0.99-1. 13//m spectral region, this implementa- 
tion of the Ackerman & Marley (2001) models predicts narrow- 
band variations of small (unndetectable) amplitude, while larger 
slope variations are predicted over 1.0-1. 08 /mi and over 1 .08— 
1.115 /mi. Large narrow -band variations, arising from the strong 
water bands, are expected redward of 1.1 2 /mi, but are diffi- 
cult to observe from the ground. In the J, H and K bands, 
Marley et al. ([2002) predict strong variations in the water bands, 
where our spectra are too noisy to draw conclusions. Variations 
similar to those of the 1.0-1.1 /urn band are expected in the 
/ band, and anti-correlated variations are expected in the K band. 
Comparison of Fig.[l0]to Fig. 7 of Ackerman & Marley (2001) 
gives a qualitative interpretation of the origin these effects. 

From these calculations we can deduce an estimate for the 
expected variations: ~3 % broad-band flux increase over 0.99- 
1.1 3 /mi and in the / band when a 5 % clear region is replaced 
by a cloudy region. According to our application of the Marley 
et al. models, narrow-band variations in the Isaac wavelength 
range and in J band, except for K i, are not predicted. In the 
H and K bands, and in Allard et al. models, the situation is 
more complex. For 2MASS J1534-2952AB, the observing con- 
ditions (high winds) and flux losses due to the slit require to 
normalise of the spectra before searching for variations (see 
Section[3]). A similar normalisation is performed on all observed 
spectra (for homogeneity) and the synthetic spectra, so that we 
are only sensitive to variations of the slope. 

We therefore calculate the variations of a spectral index, 

_ 1.00- 1.02 /mi 

S VT T — ■. 

1.06 - 1. 08 /mi 

for the data and synthetic spectra, after normalisation. The spec- 
tral index is chosen in order to maximize the predicted variabil- 
ity signal, as illustrated by the Fig.Qj] and minimize the photo- 
noise. However, small changes in the definition of the index 
(shift by +5 nm) have negligible effects on the final results, up 
to 15 % on the observed dispersion of the index (see Fig [12]) and 
on the variability upper-limits (see Table[4]l we derive. As dis- 
cussed earlier, reliable errors are necessary to properly calculate 
the upper-limits. Here we assumed that the noise from neigh- 
bouring pixels are not correlated, as is expected for photo-noise, 
and we simply square-add errors of individual pixel fluxes. 

This also applies to the IRTF observations. As a result, the 
variations we are really sensitive to, for the same cloud-cover 
variations, are now small in the J band, except at the K i dou- 
blet feature. In the H band, the main variations are due to water 
and methane absorptions at the edges of the band, mostly af- 
fecting the broad shape of the spectrum over that wavelength 
range. In the K band, variations are detected only for the coolest 
temperature (1 100K) at the red end of the band, and are due 
to variable methane absorption. Those expectations are however 
strongly model-dependent. 

For the IRTF data we use the spectral index: 

_ 1.46- 1.50/mi 
ilRTF ~ 1.58 - 1.62 /mi' 

which again maximizes our sensitivity while avoiding the main 
telluric water bands. 
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Fig. 10. Crude model predictions for the IRTF wavelength range. We used the f se j = 3 and clear models (Marley et al. 120021 in 
thick lines) and Cond, and Settl05 and Cond®2 models (Allard et al, priv.com., in thin lines). The model spectrum resolutions 
are degraded to roughly match the SpeX actual resolution. See legend of FigQT|for details. Differences between the models in the 
water bands are likely due to the different assumptions regarding the cloud opacity structure. Horizontal dashes lines indicate the 
wavelength range and effect of the normalisation in the J, H and K bands. 



Table 4. Dispersions of VLT and IRTF spectral indices, upper- 
limits (at 99 % C.L.) calculated from the observed errors, and ex- 
pected variations due to cloud cover change assuming our sim- 
ple model and a cloudy area increase of 5 % of the dwarf sur- 
face. The models are characterized by their effective tempera- 
ture and the cloud coverage, and are based on the Ackerman & 
Marley ( 120011) grid of models. 



Target 



cr(s) Upper-lim. Expected 



Model 



VLT 

2MASS J0825 
SDSSJ1254 



2MASS J 1225 
2MASSJ1534 
SDSSJ1624 



1.5% 
0.7% 



0.7% 
1.2% 
1.3% 



2.2% 
1.1 % 



2.2% 
1.6% 
2.7% 



2.1 < 
2.2' 
2.9' 
1.9' 
1.7' 
1.7' 
1.7' 



1 600 K, 95 % 
1 400 K, 50 % 
1 400 K, 95 % 
1400K, 5% 
1100K, 0% 
1100K, 0% 
1100K, 0% 



IRTF 

2MASS J0825 
SDSS 11254 



1.1% 
3.3% 



1 % 

8% 



2.6' 
4.1< 



1600K, 95% 
1 400 K, 50 % 



6.2. Interpretation 

The rotational velocity of 2MASS J0825+21 15 was measured 
by Bailer- Jones ( 120041 1 to be 16.9^gkm/s and corresponds to a 
rotational period between 2.48 and 10.75 h (with 90% of con- 
fidence levels, which include the rotation axis orientation un- 
certainty and conservative measurement errors). Hence the ob- 
serving time of each night only samples part of one period, 
while the 2003, April 16 and 17 VLT data may have been 
taken one to nine periods apart. The IRTF data, collected over 



50 hours, also cover several rotational periods. Zapatero Osorio 
et al. ( 120061 1 similarly measured vsim = 27.3 ± 2.5km/s and 
38.5 + 2.0km/s for SDSS J1254-0122 and SDSS J1624+0029 
respectively, corresponding to a period range about twice as 
short as that of 2MASS J0825+21 15. Therefore, each night 
of observations most probably covers half a rotational pe- 
riod or more of SDSS J1254-0122 and SDSS J1624+0029. 
2MASSJ1225-2739AB and 2MASS J1534-2952AB have no 
measured rotational velocity. However, brown dwarf period es- 
timates are all in the range of 1.5 to 13 hours (Basri et al. 2000; 
Moha nty & Basri [20031 Bailer- Jones [20041 Zapatero Osorio 
et al. 120061 1, so rotational periods for our targets are likely to 
fall in this range of values. If the brown dwarf's surface is a 
patchwork of clouds and cloud holes, the cloud cover ratio might 
change with time as the brown dwarf rotates. Alternatively, dy- 
namical processes could alter the ratio of cloudy to clear regions. 
In the latter case, little is known about the variability time scale 
to expect. 

We do not see the broad-band variations that would be ex- 
pected if there are large fractional changes in global cloud 
cover, using both the Ackerman & Marley ( 120011 1 and Allard 
et al. (120031 2006) models for the 1800-1 100 K effective tem- 
perature range. We show our variability upper limits, for broad 
and narrow bands, in Fig.[T3land[T4l 

We can use our spectral indices to draw quantitative conclu- 
sions. In TableH] we report the observed dispersions and 99 %- 
confidence-level upper-limits, as well as the variations expected 
based our application of the the Ackerman & Marley (2001) 
models, for a change of 5 % in the cloud cover, for the effective 
temperatures of interest. Thus, with upper-limits on the broad- 
band variability over 0.99-1. 13pm of 1.1-2.7%, depending on 
the targets, we can exclude global changes larger than ~5-8 % in 
the ratio of cloudy to clear regions, during the time we observed 
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the brown dwarfs for four of our targets. The constraint based on 
the SpeX data for 2MASS J0825+21 15 in the H band is 4 %. We 
cannot exclude the possibility that many changing regions in the 
cloud cover, each of smaller typical scale, affect a large surface 
of the brown dwarfs. 

Some of the spectroscopic time series do exhibit high- 
spectral frequency variability, which we report in Table[5] 
Generally, the amplitudes of the variations are small com- 
pared to the noise, so that we cannot confidently claim that 
these variations are real and intrinsic to the source. The low- 
amplitude variations observed at 0.997 /im in the spectra of 
SDSSJ1254-0122 and 2MASS J0825+21 15, and at 1.033/xm 
for SDSS J1624+0029, fall within an FeH band. However we de- 
tect no variations at the 0.99 /mi band head as may be expected if 
FeH absorption variations are responsible for the 0.997 /mi vari- 
ability. A search for variability at other wavelengths at which 
the object is bright and FeH is a major absorber may be indi- 
cated, such as the complex of FeH features in H band (Cushing 
et al. 120031 1. In the optical, such variations might have been de- 
tected around the CrH and FeH bandheads between 8600 and 
8700 A (Kirkpatrick et al. |2001| l. Since different bands probe 
different depths in the atmosphere, variations between the am- 
plitude of variability in the FeH bands could reveal the un- 
derlying source of the variability we observe, if it is indeed 
real. It should be noted that some low-amplitude variations 
are detected in the spectra of several targets: Around 1.03 /mi 
for SDSS J1624+0029 and SDSS J1254-0122, or 1.046 fim for 
2MASSJ1534-2952AB and SDSS J1254-0122. As we used 
different flat-field and the calibration slightly varies, we do 
not expect an instrumental effect. However we cannot relate 
those features with a known absorber. For all targets except 
SDSS J1254-0122, we do not consider that these variations im- 
pair our conclusions regarding the maximal cloud size. 

6.3. The case of SDSS J1 254-01 22 

The case of SDSS J1254-0122 is different, as in both data set, 
its spectrum exhibits numerous significant variations. A sub- 



set is common to both instruments, and a significant proportion 
of variable lines show correlated variations. For this object we 
therefore conclude that the detected variations are possibly real. 
Some, but not all, of the KI doublet lines are affected, as well as 
other lines which our simple models suggest could be variable. 
However, other features in the continuum are also seen to vary, 
which the simple model does not predict. We note that many of 
the variable features are found in region of the spectrum affected 
by telluric water band absorption. Although they sometimes cor- 
relate with each other, there is no systematic correlation, and 
other variable features are found where there are little or no tel- 
luric absorption. We also point out that some highly variable fea- 
tures detected in the SpeX data fall close the order edges or in 
regions of small SNRs, where our error determination may be- 
come inaccurate. 

Among our targets, SDSS J1254-0122 represents the proto- 
typical case of a transition object. It is therefore worth noting that 
it shows higher levels of variability than our other four targets, 
and is possibly the only one that we could classify as variable. 
This is somewhat counterbalanced by the facts that we obtained 
more good-quality data on that object, and that it is the bright- 
est (along with 2MASS J1534-2952AB). So that our sensitivity 
is higher. Finally, the small size of our sample precludes any 
conclusion regarding the variability frequency among L/T tran- 
sition objects, in comparison with earlier and later brown dwarfs. 
Nevertheless we suggest additional monitoring of this object. 

7. Conclusions 

We have obtained high SNR, low resolution spectroscopic time 
series over 0.99-1 . 1 3 //m for five brown dwarfs of spectral types 
L6 to T6, sampling the L/T transition, and for a sub-set of two 
brown dwarfs of types L6 and T2 in the J, H and K bands. We 
cautiously report some high-frequency (narrow-band) variations, 
which we generally cannot tie to specific absorption lines. An ex- 
ception is FeH in the spectra of three targets. SDSS J1254-0122 
shows numerous variable features. All require confirmation by 
additional observations. 
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Table 5. Partial list of features which we report to vary with a confidence level better than 99 %. We only report the most significant, 
or most robust varying features. We indicate in the third column additional features whose variations are correlated in time (at 99.9 % 
C.L.). 



IdTgcl 


iciescope 


waveiengin 
(fim) 


width 
(nm) 


pcdK- lO-pcd.K 


COlTcldlcll lcdllllcS 


2MASSJ0825+2115 (L6) 


VLT 


0.996 


1 


n% 






VLT 


1.008 


5 


5% 






VLT 


1.065 


5 


14% 




2MASSJ1225-2739AB (T6) 


VLT 


1.104 


1.5 


6% 




SDSSJ1254-0122 (T2) 


VLT 


0.997 


1.2 


6% 






VLT 


1.031 


3 


4% 


1.053, 1.103, 1.115 




VLT 


1.046 


4 


6% 






VLT 


1.12 


20 


10% 






IRTF 


1.127 


4 


60% 


1.472, 1.764 




IRTF 


1.471 


3 


45% 


2.324 




IRTF 


1.580 


5 


14% 


2.062 




IRTF 


1.617 


3 


25% 


2.247, 2.281 




IRTF 


2.062 


4 


22% 


2.415 




IRTF 


2.246 


5 


45% 


2.410 




IRTF 


2.281 


5 


50% 


2.324 


2MASSJ1534-2952AB (T5) 


VLT 


1.046 


4 


4% 






VLT 


1.068 


6 


7% 




SDSSJ1624+0029 (T6) 


VLT 


1.033 


6 


7% 





We place constraints on broad-band spectroscopic vari- 
ability at the levels of 2% (2MASS J1534-2952AB) to 
3% (2MASSJ0825+2115, 2MASS J1225-2739AB and 
SDSS J1624+0029), with a 99 % confidence level. 

When comparing to the variations expected from crude at- 
mospheric model interpolations between cloudy, L-type and 
clear, T-type atmospheres, we find that over the course of our 
observations, no significant variations (larger than ~5-8 % of the 
brown dwarf's disk) in the cloud coverage occurred in our sam- 
ple. Our data do not rule out smaller-size heterogeneity on the 
brown dwarf surface, nor does our analysis try to constrain other 
variability mechanisms. 
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1 1.2 1.4 1.6 1.8 2 2.2 2.4 



A (/iim) 

Fig. 1. Top: Isaac spectra of April 16 (April 17 for 
SDSS J1624+0029) for our targets, ordered by spectra type. 
The spectra are wavelength and flux calibrated using an AO 
star, smoothed over 1". The spectra are offset for clarity, as in- 
dicated in the figure. The smoothed IRTF SDSS J1254-0122 
stacked spectrum is superimposed. Bottom: SpeX spectra av- 
eraged over all nights for our two IRTF targets. Note the dif- 
ferent vertical scales for 2MASS J0825+21 15 (right labels) and 
SDSS J1254-0122 (left labels). We identify the most significant 
absorption features. 
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2MASS J0825+2115 




1.04 1.06 1.08 

A (mid) 

Fig. 2. VLT spectroscopic time series for 2MASS J0825+21 15. 
The stacked spectra are relative to the average of each night, 
and offset for clarity. UT is the mean time since April 16, 2003, 
UT=0:00. The 00:26:53 and 24:26:54 spectra have 20-min ex- 
posure time, while all other spectra have 15-min exposure time. 
The spectra are normalised over 1.014-1.054yum. The brown 
dwarf flux variations relative to the comparison star are shown 
as a black (noisy) line. The comparison star flux variations are 
in grey (almost fiat) line. The variations are relative to the inte- 
grated flux of the brown dwarf and the comparison star over that 
wavelength range. The dashed grey lines mark the features that 
we report to vary at medium resolution while the solid grey lines 
refer to the lower resolution analysis. Top panel: the sky raw 
spectra, for both observing nights. Bottom panel: the target raw 
spectra, for both observing nights. The thin-lined bluer spectra 
are those of the comparison star. None are flux-calibrated. All 
panels: all fluxes are in the same arbitrary unit, except the com- 
parison star. 
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Fig. 6. Correlation in variations for pairs of features in the 
spectrum of SDSS J1254-0122, in the VLT spectra (top) and 
IRTF (bottom). The central wavelength and the width of the nar- 
row bands are indicated on the axis. Some statistical parame- 
ters are given in the figures (particularly P, the probability that 
the variations are uncorrected). For each telescope, those pairs 
are among the most correlated ones. Ignoring any single point 
does not remove the correlation, except the bottom left one in 
the IRTF figure. 
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2MASS J1534-2952 




A Cum) 

Fig. 7. Spectroscopic time series for 2MASS J1534-2952AB. 
The 15-min spectra are relative to the average of each night. 
The flux variations of the brown dwarf (solid, black line), of the 
second comparison star (dashed, grey line), both relative to the 
bright comparison star, and the variations of the bright compari- 
son are shown (solid, grey line). All spectra are normalised over 
the 0.997-1. 054/mi. " See legend of Fig. H 

for more details. 



2MASS J1225-2739 
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Fig. 8. VLT spectroscopic time series for 
2MASS J1225-2739AB. The stacked spectra are relative 
to the average of each night, and offset for clarity. UT is 
the mean time since April 16, 2003, UT=0:00. All spectra 
have 15-min exposure time. The spectra are normalised over 
1. 002-1. 076yum. See legend of Fig.[2]for details. 
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Fig. 9. Spectroscopic time series for SDSS J1624+0029. The 15- 
min spectra are relative to the average of the night, and offset 
for clarity. All spectra are normalised over the 1. 044-1. 072/im 
region. See legend of Fig.fjfor more details. 
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Fig. 11. Crude model predictions for the VLT wavelength range, 
for three effective temperatures: 1600, 1400 and 1 1 00 K from top 
to bottom, for a gravity of 10 3 m/s 2 . We used the f sec i = 3 and 
clear models (Marley et al. 2002 ). The model spectrum resolu- 
tions are degraded to roughly match the Isaac actual resolution. 
We show the ratio of the model spectra with clouds covering 5 % 
more of their disk compared to those expected for each effective 
temperature (thick lines) and other cloudy- to total area percent- 
ages (as indicated in the graph). The amplitude of the computed 
variation depends linearly on the assumed fractional change of 
cloud coverage. 
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Fig. 12. Spectral index from the VLT (top) and IRTF (bottom) 
data. For clarity, the indices are normalised to 1 on average, inde- 
pendently for each target. (The normalisation factor is primarily 
a function of the spectral type.) Target names are indicated. 
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Fig. 13. Variation upper limits for the VLT sample, rela- 
tive to the night average (99 %-confidence-level), after smooth- 
ing over 1 1 pixels or 3 nm (thin lines), or over 5 1 pixels 
or 15 nm (thick lines). For SDSS J1254-0122 (almost iden- 
tical) independent results from both nights are shown; for 
2MASS J0825+2115, significant differences appear, due to the 
different weather conditions, and we represent the upper-limits 
mean of the two nights. Large upper-limits (redward of 1.1 //m 
for 2MASSJ1225-2739AB and SDSS J1624+0029) are trun- 
cated for clarity. 



